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ABSTRACT 

Galaxy-galaxy lensing uses the weak distortion of background sources to measure the mean excess surface 
density profile, AE(r), around a sample of foreground lensing galaxies. We develop a method for combining 
AE(r) with the galaxy-galaxy correlation function £ gg (r) to constrain the matter density parameter fl m and 
the matter fluctuation amplitude erg, going beyond the widely used linear biasing model to reach the level of 
accuracy demanded by current and future measurements. We adopt the halo occupation distribution (HOD) 
framework, and we test its applicability to this problem by examining the effects of replacing satellite galax- 
ies in the halos of a smoothed particle hydrodynamics (SPH) simulation with randomly selected dark matter 
particles from the same halos. The difference between dark matter and satellite galaxy radial profiles has a 
~ 10% effect on AE(r) at r < 1 /i _1 Mpc. However, if radial profiles are matched, then the remaining impact 
of individual sub-halos around satellite galaxies and environmental dependence of the HOD at fixed halo mass 
is < 5% in AE(r) for 0.1 < r < 15 /r'Mpc. We develop an analytic approximation to AE(r) for a specified 
cosmological model and galaxy HOD, improving on previous work with more accurate treatments of halo bias 
and halo exclusion. Tests against a suite of populated TV-body simulations show that the analytic approximation 
is accurate to a few percent or better over the range 0.1 < r < 20 /r'Mpc. We use the analytic model to in- 
vestigate the dependence of AE(r) and the galaxy-matter correlation function £ gm (r) on tt m and a%, once HOD 
parameters for a given cosmological model are pinned down by matching £ gg (r). The linear bias prediction that 
£gm( r )/£gg( r ) = constant is accurate for r > 2 /z~'Mpc, but it fails at the 30-50% level on smaller scales. The 
scaling of AE(r) with cosmological parameters, which we model as AE(r) oc fiJ^Vg , approaches the linear 
bias expectation a = j3 = 1 at r > 10 /r'Mpc, but a and (3 vary from 0.8 to 1.6 at smaller r. We calculate a 
fiducial AE(r) and scaling indices a(r) and (3{r) for galaxy samples that match the observed number density 
and projected correlation function of Sloan Digital Sky Survey galaxies with M r < -20 and M r < — 21 . Galaxy- 
galaxy lensing measurements for these samples can be combined with our predictions to constrain Q, m and erg, 
taking full advantage of the high measurement precision on small and intermediate scales. 
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1. INTRODUCTION 

In the current paradigm of structure formation, galaxies 
form by the dissipative collapse of baryons in halos of cold 
dark matter (CDM). Understanding the relation between the 
galaxy and dark matter distributions is the key challenge in 
interpreting the observed clustering of galaxies. Large area 
imaging surveys have provided a new tool for untangling this 
relationship, galaxy-galaxy weak lensing, which uses the sub- 
tle distortion of background galaxy shapes to measure av- 
erage mass profiles around samples of foreground galaxies. 
The last few years have seen rapi d growth in thi s field, with 
the first tentative detections ( Braine rdet alll9 96) giving way 
to high signal-to-noi se ratio measurements over a subs tantial 
dynamic range (e.g., |FischgxeLay|2iX^ 
iHoekstra et"aUl2002b Isheldon et alJl2004 iMandelbau m et alJ 

1 Department of Astronomy, The Ohio State University, 140 West 
18th Avenue, Columbus, OH 43210; jaiyul@astronomy.ohio-state.edu, 
dhw@astronomy.ohio-state.edu 

- Kavli Institute for Cosmological Physics, University of Chicago, 
Chicago, IL 60637; tinker@cfcp.uchicago.edu 

3 School of Natural Sciences, Institute for Advanced Study, Einstein Drive, 
Princeton, NJ 08540; zhengz@ias.edu 

4 Hubble Fellow 

5 Departments of Physics and Astronomy, University of Massachusetts, 
Amherst, MA 01003; nsk@kaka.phast.umass.edu 

6 Steward Observatory, University of Arizona, Tucson, AZ 85721; 
rdave @ as . arizona. edu 



2006). 

In a cosmological context, the strength of the galaxy- 
galaxy lensing signal for a given galaxy sample should depend 
mainly on the mean matter density ft m and the amplitude of 
dark matter fluctuations erg, since increasing either parameter 
enhances the average amount of dark matter around galaxies 
and thereby amplifies the lensing signal. 7 In this paper, we 
develop tools for constraining Q, m and erg with galaxy-galaxy 
lensing and galaxy clustering measurements, using halo oc- 
cupation models of galaxy bias that are applicable from the 
linear regime into the fully non-linear regime . Our approac h 
extends and complements e ar lier work by ISeliakl ( 2000), 
[ Guzik & Seliakl 120011 120021) . iTasitsiomi et alJ fc004 . and 
IMandelbaum et alJ <2005h 

Galaxy-galaxy lensing measures the profiles of mean tan- 
gential shear around galaxies. With knowledge of source and 
lens redshift distributions, this tangential shear can be con- 
verted to excess surface density, 



AE(r) = E«r)-E(r), 



(1) 



where E(< r) is the mean surface density interior to the disk 
of projected radius r and E(r) is the averaged surface density 
in a thin annulus of the same radius ( Miralda-Escude 1991a; 

7 Here ug is the rms linear theory matter fluctuation in spheres of radius 
8 /r'Mpc, with h = ffo/100 kms^'Mpc" 1 . 
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LSheldon et al] l2004). 8 The excess surface density profile is 
itself related to the galaxy-matter cross-correlation function 
£gm by 

AE(r) = Pc n m x (2) 
- / / r't gm Ur* + z 2 ) dzdr 1 - / ^ m (r,z)dz , 

/ JoJ-oo K ' J-oo 

where p c is the critical density of the universe. iTohnston et al.l 
( 120051) discuss and test methods of inverting A£(r) to obtain 
the three-dimensional £ gm (r). Here we treat AS(r) as the pri- 
mary observable and concentrate on predicting it directly. 

On large scales, where matter fluctuations are linear, the 
relation between the matter auto-correlation function £ mm , the 
galaxy-matter cross-correlation function £ gm , and the galaxy 
auto-correlation function £ gg may be adequately described by 
the linear bias model, 

£gg = b (mm, (3) 
£gm = ^£mmi (4) 

where the lin ear bias factor b is the same in both equations 
jKaiseJll984l) . Thus, measurements of £ gg and AS oc £ gm f2 m 
can be combined to yield fl m /b. Since the amplitude of galaxy 
fluctuations is crg ig = /?erg in the linear bias model, this method 
in turn constrains the product <7gf2 m . Redshift-space distor- 
tions of the galaxy power spectrum and the abundance of rich 
galaxy clusters as a function of mass both depen d on a param- 
eter combination that is approximately crgSl^ 6 (Kaisen fl987t 
IWhite et all 19931) . so the combination of galaxy-galaxy lens- 
ing with either of these measurements can break the degener- 
acy between erg and ft m . However, the linear bias approxima- 
tion may break down on the scales r < several /r'Mpc where 
AS(r) is measured with high precision. Moreover, if the rela- 
tion between galaxy and matter density contrasts is linear but 
stochastic, then the linear bias factor b in equation @ should 
be replaced by br„ m , where r„ m is the galaxy-m atter cross- 
correlation coefficient (lPenlll998t lDekel & Lahav 1999), and 
the constrained combination becomes a%Q, m r Spl even in the 
linear regime. The addition of r sm as a free parameter re- 
duces the cosmological constraining power of the AE and £ gg 
combination, and restoring it requires an independent deter- 
mination of £ mm . Co smic shear measurements can provide 
such a determination ( Blandfo rd et aljH99ll Miralda-Escude 
1991bi lKaiserll9 92). but these measurements are challenging 
a nd suffer larger systema tic errors than galaxy-galaxy lensing. 

iHoekstra et all 120011 [2002) used imaging and photomet- 
ric redshift data from the Red-Sequence Cluster Survey 
( Glad ders & Yedl2001l) to measure aperture fluctuations pro- 
portional to £ gg , £ gm , and £ mm . They provided tentative evi- 
dence that b and r are each individually scale-dependent but 
that the ratio b/r gm is approximately constant, with b/r gm ~ 1 
for tt m = 0.3. U sing the Sloan Digital Sky Survey (SDSS, 
lYork et aTll2000t) . which provides spectroscopic redshifts of 
lens galaxies and image shapes and pho tometric redshifts 
of source galaxies, She ldon et al.l ll2004l) detected galaxy- 
galaxy lensing and measured the galax y-matter corre lation 
fun ction from 0.025 to 10 /r'Mpc (see Fischer et all 120001 
and McKav et al. 2001 for earlier SDSS measurements, and 
Mandelba um et al.l 12006 for more recent measurements at 
r < 2 /r'Mpc). They found that the galaxy-galaxy cor- 
relation function and the galaxy-matter correlation function 

8 We interchangeably use r to refer to a projected (two-dimensional) or a 
three-dimensional radius. 



agree in shape, with an amplitude ratio that implies b/r gm = 
(1.3 ±0.2)(f2 m /0.27) for galaxy samples of mean luminosity 

To circumvent the limitations of the linear bias approxima- 
tion, we model galaxy-galaxy and galaxy-matter correlations 
using ha lo occupation method s , fol lowing the lead of Seliak 
(2000), iBerlind & Weinberg! (120021) . and iGuzik & Seliak 
(2002). The halo occupation distribution (HOD) provides 
a fully non-linear description of the relation between galax- 
ies and mass by specifying the probability P(N\M) that a 
halo of virial mass M contains N galaxies of a particu- 
lar class, along with any spatial or velocity biases within 
individual halos. 9 Numerous authors have used this ap- 
proach to compute analytic approximations fo r galaxy and 
dark matter clustering stat istics ("e .g.JMa & FrvK OOO: Seliak 
1 20001 IScoccimarro et all 1200 It ISeliakl 1200 it IShefh et alJ 
1200 1 at IWhitd 1200 ll see review bv ICoorav & Shefhll2002ir 
and to model observed galaxy clusteri ng ("e.g.. [j ing et al] 
1 19981. 120021 iPeacock & Smithl 120001 iKochanek & White! 
[ 20011 iBullock et alJ 120021 iM agliocchetti & Porcianil |2003[ 
l Yang et all 120031 iPorciani et alJ 120041: IZehavi et all 120041 
2005aUZhengl2004tlA bazaiia n et al.l2005tlColiister & Lahavl 
l2005tlLee et all2006UTinker et all2 006). Theoretical predic- 
tions for the HOD of different galaxy types have been calcu- 
lated using semi-analytic models, hydrodynamic simulations, 
and high resolution iV-body calculati ons that identify "galax- 
ies" with dark mat ter substructures llK auffmann et al]|1222 ; 
Benson et alJ2000[ [White et alJ2001l:lYoshikawa et alJ | 2001 ; 
Berlind et alJl2003UKravtsov et alJl2004HZentner et al.ll2005 : 
Zheng et 31120051 

Our basic approach to modeling galaxy-galax y lensing in 
the HOD fra mework is similar to that adopted by Tink er et"all 
(2005, 2006) for modeling mass-to-light ratios and redshift- 
space distortions. Since measurements of the galaxy power 
spectrum, cosmic microwave background anisotropies, and 
the Lya forest yield tight constraints on t he shape of the linear 
matter power spectrum P\\ n (k) (see, e . g..lSpe rgel et al. ]|2003l 
Tegmarket al. 2004: ICole et al.ll2005l Iseliak et alJl2005l) . we 
take this shape to be fixed and investigate the parameter space 
spanned by 0,„ and erg. For a given choice of Vl m and erg, 
we first choose HOD parameters to match observations of 
the projected galaxy corr elation function w p (r p ) (see, e.g., 
Zehav iet alJl2004l 12005 al) . We then predict the excess sur- 
face density profile AE(r) for this combination of f2 m , a%, and 
HOD. Comparison to galaxy-galaxy lensing measurements 
then determines the acceptable combinations of f2„, and erg. 
We impose the observed galaxy-galaxy correlation function 
as a constraint on the HOD, instead of taking ratios as in 
the linear bias analysis. There is no need for an unknown 
cross-correlation coefficient r gm because any "stochasticity" 
between galaxy and mass density fields is automatically incor- 
porated lm_tJie_IiOr2ralculadori i _Our str ategy co mplements 
that of IGuzik & Seli aid ( 1200 ll 120021) and Mandelbaum et al. 
(2005), who focus on constraining halo masses, halo pro- 
files, and satellite fractions using AE(r) alone, rather than 
constraining VL m and erg from the combination of AE(r) and 

CggW. 

Our eventual conclusions about the cosmological constrain- 
ing power of galaxy-galaxy lensing measurements rest on an 
analytic model for computing A£(r) given Pu n (k), Q, m , erg, 

9 Throughout this paper, we use the term "halo" to refer to a dark matter 
structure of overdensity p/p m ~ 200, in approximate dynamical equilibrium, 
which may contain a single bright galaxy or a group or a cluster of galaxies. 
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an d the galaxy HOD. Th is model is similar in spirit to that 
of iGuz ik & Seliak (2002), but it differs in many details, in 
part because we define the calculational problem in different 
terms. We test the analytic model against numerical calcula- 
tions, in which we use a specified HOD to populate the halos 
of A^-body simulations, placing "central" galaxies at the halo 
potential minimum and "satellite" galaxies at the locations of 
randomly selected dark matter particles. Both our analytic 
model and our method of populating A^-body halos ignore the 
impact of dark matter subhalos around the individual satellite 
galaxies orbiting in a larger halo. To begin, therefore, we test 
the validity of the "populated halo" approach itself, by com- 
paring A£(r) for the galaxy population of a smoothed particle 
hydrodynamics (SPH) simulation to that found by populat- 
ing the dark matter halos of this simulation with "galaxies" 
placed on randomly selected dark matter particles. This test 
shows that satellite subhalos have minimal impact on A£(r) 
and that the populated halo approach is acceptable for our 
purpose. We also show that any environmental dependence 
of halo galaxy co ntent at fixed halo mass (iGaoetal J 120051 
lHarker et alJ200 6) has little discernible impact on the galaxy- 
galaxy or galaxy-matter correlation functions in our SPH sim- 
ulations. More generally, our SPH and A^-body tests indicate 
that the analytic model should be accurate at the 5-10% level 
on scales r > 0.1 /r'Mpc. This level of accuracy is acceptable 
for present purposes, since the curre nt measurement error s are 
typically > 25% per radial bin (e.g., Sheldon et al. 2004J), but 
still higher accuracy will be needed in the long term. 

In our A^-body and analytical calculations, we use HOD pa- 
rameters for SDSS galaxy sam ples with abs olute-magnitude 
limits M r < -20 and M r < -21 JZehavi et alJl2005 b) for pur- 
poses of illustration. 10 The results presented in §|5]therefore 
yield predictions of the weak lensing signal for these galaxy 
samples as a function of Vt m and a%. The analytic model can 
be used to make predictions for other galaxy samples given 
measurements of the projected correlation function as input. 

2. SPH GALAXIES VERSUS POPULATED HALOS 

To test the validity of our A^-body method for calculating 
galaxy-galaxy lensing predictions (see § |3j, we first examine 
an SPH simulation of a ACDM (inflationary cold dark mat- 
ter with a cosmological constant) universe. This simulation 
is described in detail bv lWeinberg etalJ {2004), who, among 
other things, present predicted galaxy-matter correlations and 
compare them to recent observations. Here we want to know 
whether the individual dark matter subhalos retained by bary- 
onic galaxies in groups and clusters make an important con- 
tribution to the galaxy-galaxy lensing signal. 

In brief, the simula tio n uses Parallel Tre eSPH 
JHernquist & KatzJ Il989t iKatz. Weinberg. & Her nauist 
1996; iDave et al.l 1 1 9971) to model a 50 /r'Mpc comoving 
cube with 144 3 dark matter particles and 144 3 gas particles. 
The cosmological parameters are Vt m = 0.4, fl\ = 0.6, 
h = 0.65, n = 0.95, Q b h 2 = 0.02, and cr 8 = 0.80. The grav- 
itational forces are softened with a 10 /r'kpc (comoving) 
spline kernel. Radiative cooling leads to the formation o f 
dense baryonic clumps (IKatz et alJll992HEvrard et alJI 19941) . 
which form stars according to the algorithm described by 
IKatz. W einberg ^ Hernquisl (119961) . Galaxies are identified 
by applying the SKID (Spline- Kernel I nterpolated DEN- 
MAX; see IKatz. Weinberg. & Hernquistl 119961) algorithm 



to the population of stars and cold, dense gas particles. 
Tests with simulations of varying resolution show that the 
simulated galaxy population is complete above a baryonic 
mass (stars plus cold, dense gas) of ~64 wsph, corresponding 
to 5.4 x 10 10 M Q (3.5 x 10 1() /r'M ) for this simulation. 
The space density of galaxies above this mass threshold is 
n g = 0.02 (/i~'Mpc)~ 3 , corresponding to that of observe d 
galaxies with M r < -18.6 (L > 0.18 Z^: lBlanton et alJl2003l) . 
We use this mass-thresholded galaxy sample for the tests 
below. 

We identify dark matter halos by applying the friends-of- 
friends algorithm (FOF; Dav is et allll9 85) to the dark mat- 
ter particle distribution, with a linking length of 0.2 times 
the mean interparticle separation, or 70 /r'kpc. We associate 
each SPH galaxy with the halo containing the dark matter par- 
ticle closest to its center of mass. To create "populated halo" 
galaxy catalogs, we replace the SPH galaxies in each halo 
with an equal number of artificial galaxies positioned on dark 
matter particles. The first "central" galaxy of each occupied 
halo is placed at the location of the dark matter particle with 
lowest potential energy (computed using only halo members). 
Any additional, "satellite" galaxies are placed on randomly 
selected dark matter particles. Satellites therefore follow the 
radial profile of dark matter within each halo, while any de- 
tailed association between satellites and the centers of dark 
matter subhalos is erased. 

The left panels of FigureEJcompare the galaxy-galaxy cor- 
relation functions, galaxy-matter correlation functions, and 
excess surface density profiles of the SPH galaxies and the 
populated halo galaxy catalogs. Results for the populated ha- 
los are an average over 10 realizations of the galaxy locations, 
and error bars show the dispersion among the ten realizations. 
(Note that these do not represent the uncertainty on the mean, 
which would be a factor of three smaller, and they do not in- 
clude the uncertainty owing to the finite simulation volume, 
since we are comparing galaxy catalogs in the same volume). 
The galaxy-galaxy and galaxy-matter correlations of the two 
catalogs are very similar at r > 0.5 /r'Mpc, while at smaller 
separations the populated halo catalog has correlations that 
are stronger by up to 20%. The excess surface density profile 
is calculated by directly counting galaxy-dark matter parti- 
cle pairs in projection to compute £(< r) and S(r), not by 
integrating the three-dimensional £ gm (r). We count all pro- 
jected pairs through the 50 /r'Mpc box and average results 
from the three orthogonal projections of the box; noise from 
uncorrelated foreground and background particles cancels out 
because we have many galaxy targets. Relative to the SPH 
galaxy catalog, AS(r) for the populated halo catalog starts 
about 10% low, rises to 10% above at r ~ 0.5 /r'Mpc, then 
agrees closely beyond r ~ 1 /r'Mpc. 

The modest deviations between the SPH galaxy and popu- 
lated halo results could reflect either the impact of satellite 
subhalos or differences between the radial profiles of SPH 
satellites and dark matter. To separate the two effects, we 
adopt a different method of populating halos that ensures 
identical radial profiles, by placing each satellite at the ra- 
dial distance of the corresponding SPH satellite but choos- 
ing a random orientation for the radius vector. Results are 
shown on the right panels of Figure [2 The differences in 
£ gg (r), £ gm (r), and A£(r) are greatly reduced, demonstrating 
that they arise mainly from the different radial profiles of SPH 



Throughout the paper, we quote absolute magnitudes for h = 1; more 
generally, these thresholds correspond to M r — 5 log h. 



11 We use the implementation of SKID by J. Stadel & T. Quinn, available at 

http : / / www-hp cc .astro. Washington. edu/tools/skid .html 
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FIG. 1. — Galaxy-galaxy correlation functions (top panels), galaxy-matter correlation functions (middle panels), and AS(r) profiles (bottom panels) for the 
true galaxy population of an SPH simulation (solid lines) and the populated dark matter halos of this simulation (dotted lines, see text). Inset panels show the 
fractional difference between the SPH and populated halo results. In the left-hand panels, satellite galaxies in the populated halos are placed on randomly selected 
dark matter particles, while in the right-hand panels they are forced to follow the radial profile of satellite galaxies in the SPH simulation. 
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FIG. 2. — Comparison of galaxy-matter correlation functions for satellite 
galaxies of our standard SPH simulations (dashed lines), which uses 144 3 
particles in a 50 /r'Mpc box, to those of a higher resolution, smaller vol- 
ume simulation (128 3 particles in a 22.222 /r'Mpc box, solid lines). In 
both simulations we select halos in the mass range 6 X 10 i2 /!~'Mq < M < 
2 X lO 13 /! _1 A/0 and satellite galaxies above the 3.5 X 10 10 /i~'Mq resolu- 
tion limit of the larger volume run. Heavy lines show the full satellite-matter 
cross-correlation function, while light lines include only matter in the satel- 
lite's parent halo (the one-halo term). Error bars are computed for the smaller 
simulation via bootstrap resampling. 



satellite galaxies and dark matter; specifically, the SPH satel- 
lites are less concentrated towards the halo center than the 
dark matter. With matched radial profiles, the populated ha- 
los still have slightly larger £ gm (r) at r < 0.1 /r'Mpc, in part 
because there are usually offsets of this magnitude between 
the location of the SPH central galaxy and the position of the 
most bound dark matter particle. However, the differences in 
AS(r) are now smaller than 10% at all r. 

We conclude that it is safe to ignore the subhalos of in- 
dividual satellite galaxies when computing A£(r) for a full 
galaxy sample. Indeed, the remaining residuals in Figure ^ 
a consequence of the central galaxy offsets mentioned above, 
are opposite in sign to those expected from satellite subhalos. 
Satellites in the SPH simulation do res ide in individual dark 
matter subhalos JWeinberg et al J 120061) . but these are tidally 
truncated, and at small separations A£(r) is dominated by the 
contribution of the more numerous, central galaxies (see § 14.21 
below). The small impact of satellite subhalos on AE(r) is 
therefore u nsurprising, and was anticipated by earlier analy tic 
modeling (Guzik & Selia kl2002l lMandelbau m et al.l20 05). 

Nonetheless, one might worry that the absence of any sub- 
halo signal in FigureJ^is an artifact of our simulation's finite 
mass resolution, leading to an artificially high degree of tidal 
truncation. To test this possibility, we compare results from 
this simulation to those of a simulation of the same cosmo- 
logical model with a factor of eight higher mass resolution 
but smaller volume. This simulation uses 128 3 dark matter 
particles and 128 3 gas particles in a volume 22.222 /r'Mpc 
on a side. In each simulation, we select all halos in the mass 



range 6 x 10 12 /t'Mq < M < 2 x 10 13 r'M and measure 
the galaxy-matter correlation function for satellites in these 
halos above the baryonic mass threshold of the larger volume, 
lower resolution simulation. The 50 /i _I Mpc box contains 78 
halos and 188 satellite galaxies satisfying these cuts, while 
the 22.222 /r'Mpc box contains 12 halos and 33 satellites. 
As shown in Figure|2] the satellite galaxy-matter correlations 
are equal in the two simulations to within the statistical errors, 
which are estimated by bootstrap resampling of the galaxies 
in the smaller simulation. The average mass profiles around 
satellites are therefore robust over a factor of eight in mass 
resolution. 

Standard HOD calculations assume that the halo occupa- 
tion function P(N\M) has no direct dependence on a halo's 
larger scale environment. This assumption is motivated by 
the excursi on set derivation o f the Extended Press-Schechter 
formalism ( Bond et an fl99lh . which, in its simplest form, 
predicts that a halo's formation history is un corrected with 
its environment at fixed mass (White 1996). The corre- 
lation of galaxy properties with large scale environment 
emerges indirectly from the correlation with halo mass be- 
cause h igh mass halo s are m ore common in dense environ- 
ments. Blanto net alJ (|2Q06) showed that the observed cor- 
relation of red galaxy fraction with overdensity measured at 
6 /r'Mpc is entirely accounted for by the correlation with 
overdensity measured at the 1 /r'Mpc scale characteristic of 
individual large halos. However, while early A^-body stud- 
ies showed at most weak correlations between halo forma- 
tion time an d environment at fixed mass for halos with M > 
10 13 /t " 'M^ JLemson & Kauffmannl Il999t ISheth & Tormenl 
2004), iGao et all (|2005) have recently shown that there is a 
much stronger correlation for lower mass halos, with the older 
halos being more strongly clustered (see also She th & Tormenl 
120021 lHarker et al J 12006 discuss the potential origin of en- 
vironmental correlations in the excursion set formalism). 
iBerlind et alJ ( 120 03). examining the same SPH simulation and 
galaxy sample that we have used here, showed that the mean 
number of galaxies as a function of halo mass, (N)m, is inde- 
pendent of halo environment within the statistical uncertain- 
ties imposed by the finite simulation volume. However, in 
light of Gao et al.'s (2005) result, we have carried out an ex- 
periment to explicitly examine the possible impact of environ- 
mental dependence of P(N\M) on galaxy-galaxy and galaxy- 
matter correlations. 

In the populated halo calculations shown in Figure [2 the 
number of galaxies assigned to each halo is equal to the num- 
ber of SPH galaxies, so any environmental dependence pre- 
dicted by the SPH simulation is also built into the populated 
halo distribution. We eliminate the environmental dependence 
by shuffling the galaxy populations among halos of similar 
mass. Specifically, we sort the halos by mass and replace the 
number of galaxies A^,- in the halo of rank i with the number 
Ni+\ in halo z"+ 1, then recalculate £ gg (r), £ gm (r), and AS(r), 
averaging over ten realizations of galaxy positions within ha- 
los. We repeat the exercise with the substitutions Nj+\ — > Nj+2, 
Ni-i, and N1-2 so that we can average over four different halo 
shufflings and compute the statistical error on the mean. The 
sampling of the halo mass function is sparse for the highest 
mass halos in the simulation, so at high masses we cannot 
exchange the galaxy contents of halos without significantly 
changing P{N\M) itself. We therefore keep the galaxy popu- 
lations of the Affi x most massive halos fixed, with Affi x = 5 or 
20. For Affi x = 5, we are shuffling the contents of all halos with 
M < 9.0 x 10 13 /r'M Q , and for = 20 we are shuffling the 
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TABLE 1 

Parameters of gadget Simulations and HOD Parameters 



Model 






AfminCr'Afo) 


Mi(h- l M Q ) 


«sat 


Avir 


1 


0.10 


0.95 


3.81 x 10" 


1.22 x 10 12 


1.04 


240 


2 


0.16 


0.90 


6.27 x 10" 


1.73 x 10 13 


1.03 


210 


3 


0.30 


0.80 


1.17 x 10 12 


2.65 x 10 13 


1.04 


190 


4 


0.47 


0.69 


1.87 x 10 12 


3.49 x 10 13 


1.05 


165 


5 


0.63 


0.60 


2.41 x 10 12 


3.85 x 10 13 


1.09 


155 



NOTE. — The HOD parameters of the fiducial model (Model 3) 
are chosen to reproduce the same clustering of the SDSS galaxy sam- 
ple of M r < -20 and to match the number density of galaxies n g = 
5.74 X 10~ 3 (/r'Mpc) -3 . HOD parameters of the other models are scaled 
with f2,„ from the HOD parameters of the fiducial model and adjusted to 
match £ gg and n g . 



contents of all halos with M < 4.6 x 10 13 h~ x M & . 

Figure|3]plots the fractional difference in £ gg (r), £ gm (r), and 
AS(r) between the shuffled halo realizations and the original 
populated halos. We use the populated halos as the compari- 
son standard rather than the SPH galaxies so that we can iso- 
late the impact of environmental dependence of P(N\M). Er- 
ror bars show the error on the mean from the four shufflings, 
but recall that we have only one realization of the original 
populated halos. For A^ x = 5, there is a 5% increase on £ gg (V) 
at r < 0.5 /i _1 Mpc. However, these scales lie in the 1-halo 
regime where environmental variation of P(N\M) should have 
no impact at all, so the increase is probably a statistical fluc- 
tuation that reflects the particular sizes and concentrations of 
the halos present in the simulation. It is only slightly larger 
than the la error bars, and the errors from point to point are 
highly correlated. For Nsx = 20, the changes in £ gg (r) are less 
than 3% over the range 0.08 /z _1 Mpc < r < 3 /r'Mpc. The 
three points at r > 5 /r'Mpc are depressed by ~ 5% on aver- 
age, which suggests that shuffling may slightly lower the large 
scale galaxy bias factor, but the statistical significance of this 
depression is difficult to assess with a single 50 /r'Mpc sim- 
ulation. 

Shuffling changes £ gm (r) by less than 5%, usually much 
less, except for the largest scale point with Nf\ x = 5. Most sig- 
nificantly for our present purposes, the changes in AE(r) are 
at most ~ 2% for A^ x = 20 at all scales, and only slightly larger 
for /Vfix = 5. We conclude that ignoring any possible envi- 
ronmental dependence of P(N\M) has minimal impact on the 
calculation of galaxy-galaxy lensing observables for a given 
cosmology and HOD, for a galaxy sample defined by a thresh- 
old in baryonic mass. There could be a few percent effect 
on the large scale bias of the galaxy-galaxy correlation func- 
tion, which might lead to small errors in inferring the HOD 
from observations of £ gg (V). Assessing the impor tance of this 
effect will require larger simulations. ICroton. Gao & White! 
(2006) have carried out a similar shuffling experiment for 
semi-analytic galaxy populati ons in the 5 00 /r'Mpc Millen- 
nium Run simulation (Croton et al. 2006), and they find few 
percent changes in large scale bias for galaxy samples defined 
by thresholds in mass or absolute magnitude (it was the hear- 
ing about their shuffling experiment that inspired us to carry 
out our own). 

3. A'-BODY SIMULATIONS 

To help us develop and test our analytic model, we have 
carried out five A^-body simulations of a ACDM universe us- 
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FIG. 3. — Possible impact of environmental variation of the HOD on £g g (r), 
£„ m (r), and AE(r). We shuffle the occupation numbers of halos of similar 
mass, leaving the populations of the five (filled circles) or 20 (open squares) 
most massive halos unchanged. Plots show the fractional difference between 
the shuffled halo results and the original results. Error bars in the points 
show the uncertainty in the mean calculated from four different shufflings 
(see text). 
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FIG. 4. — Large panels show N-body results for § gg (r), ^ gm (r), and AS(r) for the five cosmological parameter combinations indicated in the legend and 
detailed in TablelTl Attached bottom panels show the fractional difference between the analytic model calculations and the simulation results. Error bars 
represent fractional statistical uncertainty on the W-body results for the central model (f2,„ = 0.3 and ag = 0.8, solid lines), computed from the error on the mean 
of the five simulations. 



ing gadget (iSpringel et alJ l200ll) . Each simulation begins 
at expansion factor a = 0.01 with a scale-invariant (n = 1) 
fluctuation spectrum m odulated by the transfer function of 
Efstat hiou et al.l(ll992l) with shape parameter V = 0.2. Our an- 
alytic mod el calculations i n §E]use the cmbfast transfer func- 
tion (|Seliak & Zaldarriagal ll 9961) . which re presents cosmo- 
logica l predictions more accurately, but the Efstathio u et al.l 
( 1992) representation should be adequate for calibrating and 
testing the analytic model itself. The simulations end at 
a = 1.0, when il m =0.1, Ha = 0.9, and the linear theory nor- 
malization of the power spectrum is as = 0.95. We use ear- 
lier outputs from the same simulations to represent models 
with the cosmological parameter combinations listed in Ta- 
ble □ (fi m , ct 8 )=(0.16, 0.90), (0.30, 0.80), (0.48, 0.69), and 
(0.63, 0.60). Since we are adopting a fixed, observationally 
motivated form of the power spectrum instead of changing its 
shape with f2 m , this procedure is exact. We would obtain the 
same results if we ran a separate simulation for each model but 
started it at expansion factor a = 0.01 /a out , where a out = 0.84, 
0.64, 0.49, 0.40 for the four (£l m , q- s ) co mbinations. We re- 
fer the reader to iTinker et al J ]2005l 120061) for the simulation 



details. 

Our simulations use 360 3 particles to model a volume 
253 /z _I Mpc (comoving) on a side. The dark matter particle 
mass is 9.6 x 10 10 n m /z _1 M Q . We choose the mass resolution 
so that the lowest mass halos that host galaxies with M r < -20, 
according to our HOD fits (see below), contain at least 32 par- 
ticles. The gravitational force resolution is e = 70 /r'kpc (this 
is the approximate Plummer-equivalent value). The five sim- 
ulations are identical except for the random number seed used 
to generate the initial conditions. We identify dark matter ha- 
los using FOF with a linking length equal to 0.2 times the 
mean interparticle separation, or 140 /r'kpc, and set the halo 
mass equal to the total mass of the linked particles. 

We populate the A^-body halos with galaxies using HOD pa- 
rameters that are designed to reproduce the mean space den- 
sity and projected correlation function of SDSS galaxies with 
M r < -20, as measured bv lZehavi et alJ(l2005bJ). The ad opted 
form o f the HO D is mot ivated by the results o flKravtsov et alJ 
(2004) and lZhenget al.l fc005).' Halos below some minimum 
mass M m ; n are devoid of galaxies. All halos above M m ; n have a 
central galaxy, which is placed at the position of the dark mat- 
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ter particle with the lowest potential energy in each halo. The 
number of satellite galaxies is drawn from a Poisson distri- 
bution with mean (M /M\) a ™ 1 . Each satellite galaxy is placed 
on a randomly selected dark matter particle from the halo. 
Table [0 lists the values of M m i n , M\, and a sat for our five 
(f2 m , eg) combina tions. Further detai ls of the fitting proce- 
dure are given bv iTinker et al.l ( 1200 6). The specifics of the 
parametrization and details of the fitting method are not im- 
portant to our purposes here, since we will test the analytic 
model predictions using the same HOD parameters applied to 
the simulations. However, these parameter choices ensure a 
galaxy population with realistic clustering properties. 

Figure |4] shows £ gg (r), CgmW) an d AE(r) for the five N- 
body models. The five galaxy-galaxy correlation functions 
are nearly identical by construction, though with the HOD 
parameters at our disposal it is not possible to exactly match 
the observed correlation function over our full range of a%. 
The galaxy-matter correlation function is higher for the more 
strongly clustered, high a% models, as expected. However, 
since AS(r) scales (approximately) with Vt m a%, and VL„, falls 
faster than eg grows in our simulation outputs, the order of 
models is reversed on the A£(r) panel. We discuss the com- 
parison between the A^-body and analytic model results in the 
following section. 

4. ANALYTIC MODELING OF GALAXY-MATTER 
CLUSTERING 

4. 1 . Formulation and Tests 

Our analytic method of calculating £ gm (r) for a given 
cosmo logy and HOD is based on the methods that IZheng l 
(2004) and Tinker et al. (2005a, see Appendix B) 
used to calculate the galaxy-galaxy correlation function. 
These methods are based, in t ur n, on idea s introd u ced by 
Schen- er & Bertschingerl (119911). iMa&Frvl d2000l). iSeliakl 

(2000) , P eacock & Smith! ( 120001) and |Scoccimarro et all 

(2001) . We present a full technical description of our £ gm 
calculation here but refer the reader to these earlier works 
for more general d iscussion. Our gala xy-galaxy correlation 
calculations follow Tinker et al. (2005), with ellipsoidal halo 
exclusion. 

Contributions to £ gm can come from galaxy-matter pairs 12 
residing in a single halo or in two distinct halos. We separate 
these two contributions as 



i + w) = [i + Cm] + Cw] ' 



(5) 



noting that it is pair counts (proportional to 1 + £ gm ) that add 
rather than the correlations £ gm themselves. The one-halo con- 
tribution is 



i 



47r r 2 n„ 



" dn , , M 
dM—(N) M 



1 



A/„, 



dM 



-F' 



2R V 



£6) 



where dn/dM is the halo mass function (| Press & Schechterl 
H974t lSheth&TormenlfT999t Uenkins et al.l 120011) . (N)Z~is 
the mean number of galaxies in halos of mass M, p m is the 
mean mass density, and F(r/2R v i r ) is the average fraction of 
galaxy-matter pairs in halos of ma ss M and virial radius R v j T 
that have separation less than r (iBerlind & Weinberg 120021 
F'(x) is simply the derivative of F(x) with respect to its argu- 
ment). We define such that the mean density within /? v „ 
is A v ; r/ o m , and unless otherwise stated we assume A v i r = 200. 
We further split the one-halo term by discriminating central 

12 By which we mean pairs of galaxies and dark matter particles. 



and satellite galaxies (see, e .g., IBerlind & Weinberg 120021 
lYan g et a l.l2003l7.hen g |2004l). 

(N) M F'(.x) = (N cen ) M F' cen (x) + {N^) M F' s Jx). (7) 

Pairs involving a central galaxy simply follow the radial mass 
profile p(r), so F' cen (x) oc p(r)r 2 . The distribution Fj at (x) of 
satellite galaxy-matter pairs is the convolution of the galaxy 
and matter profiles. We assum e a spherical NFW profile 
llNavarro. Frenk & White! 1 19971) . truncated at R vil , for both 
dark matter and satellite galaxies. We compute the dark 
matte r concentration parameter c,j m using the relation of 
iBullock et alJ (120011) . We allow the galaxy concentration to 
be different, c ga i = a c Cd m , but adopt a c = 1 as our standard 
assumption. 

On scales much larger than the virial diameter of the largest 
halo, the galaxy-matter correlation function is equal to £ mm (r) 
multiplied by a galaxy bias factor 



baa = ■ 



' dn 

dM—(N) M b h (M), 
dM 



(8) 



where bh(M) is the bias factor of halos of mass M. How- 
ever, an accurate calculation on intermediate scales must ac- 
count for the finite extent of halos, for the scale dependence 
of b^iM), and for halo exclusion — two spherical halos can- 
not be separated by less than the sum of their virial radii. It 
is convenient to do the calculation in Fourier space, where the 
convolutions of halo profiles become multiplications of their 
Fourier transforms. Our complete series of expressions for the 
two-halo contribution to ^gmW is 



i+C'w 



i 



where 

C'W = j^ J J kk 1 Pitkin 
is the Fourier transform of 



(9) 



(10) 



1 Z" 00 ✓/ 

P^(k\r) = P m (k)- / dMi — {N) M MMi\r)y g (k,Mx) 
n' e J dM\ 

b h (M 2 \r) y m (k,M2)Pn (x\M u M 2 ), (11) 



s jo 

f°° dn M 2 

x / dM 2 

J dM 2 p m 

where y g (k,M) and y m {k,M) are the normalized Fourier coun- 
terparts of the galaxy and the matter profiles, and 



dM, 



dn 
d~M\ 



(AO 



A/i 



J dn , , 

dM 2 ——(N) M2 p ao (xM u M 2 ). 
dM 2 



(12) 

In these expressions, p ao {x\M\ : M 2 ) represents the probabil- 
ity that two halos of mass M\ and M 2 with scaled separa- 
tion x = r/(/?vir,i +^vn,2) do not overlap. For spherical ha- 
los, Pn n(x) would be a step function at x = 1, but Tink er et"all 
(2005) found that an accurate separation of the 1-halo and 
2-halo contributions to £ gg (r) requires accounting for the non- 
spherical shapes of halos identified by the FOF algorithm. We 
adopt their expression, based on a fit to Monte Carlo realiza- 
tions of ellipsoidal halo pairs with a reasonable distribution of 
axis ratios: p ao (x) = forx < 0.8, p ao (.x) = 1 for jc > 1.09, and 
/7 no (jc) = (3y 2 -2y 3 ) with y = (x-0.8)/0.29 for 0.8 <x< 1.09. 
The restricted number density n' is the mean space density of 
galaxies residing in allowed (i.e., non-overlapping) halo pairs 
at separation r. Since p no (x\Mi,M 2 ) and the halo bias factors 
bh(M\r) depend on r, one must evaluate equations (I9l-( I12> 
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separately for each value of r where one wants to know £ grn (r). 
The double integrals in equations (II li and d!2t are non- 
separable because of the Mi -dependence of p no (x). In prin- 
ciple, one should separately compute equation il l\ for cen- 
tral and satellite galaxies and sum the results, since y g (k,M) is 
different in the two cases, but we have tested and found that 
ignoring this subtlety has negligible effect. 
For scale-dependent halo bias factors, we adopt the expres- 



1.49 



b 2 h (M\r) = b 2 (M)x 



[l + 1.17g mm (r)] 

n 2.09 



(13) 



[l+0.696 nm (r)r 

from iTinker et all (12001 . We also use Tinker et al.'s 
(2005a; Appendix A) expressions for the asympto tic bias fac- 
tors fras vmCM). These follow the formulation of ISheth et al.l 
(2001b), but with different parameter values that yield a sub- 
stantially better fit to the simulations. For "concordance" 
cosmologi cal parameters, these bias factors are similar to 
those of Seli ak & Warren (2004), but they are more accu- 
rate for models with different matter power spectra. We use 
Smith et al.'s (2003) approximation for the non-linear power 
spectrum P mm (k) and correlation function £ mm (r) in equa- 
tions (1111 and (II 31 . One could in principle use P\i n (k) in- 
stead of P mm (k); this would require a different (though still 
scale-dependent) expression for bh(M\r) with a separate N- 
body calibration. 

The calculation as we have described it is a straightfor- 
ward generalizatio n of the £ gg (r) calculation presented by 
ITinker et alJ (|2005), who tested its accuracy over the range 
a% = 0.6-0.95 for both the T = 0.2 simulations used here and 
for a similar set with T = 0.12. However, a significant techni- 
cal, and to some degree conceptual issue arises with the evalu- 
ation of the second integral in equation (II It . Since we assign 
all galaxy-matter pairs to either the 1-halo or 2-halo terms, 
we implicitly assume that all dark matter is in halos of some 
mass, and thus 



dn 

dM 

dM 



M = p„ 



(14) 



More importantly for present purposes, the distribution of 
matter is by definition unbiased with respect to itself, and 
therefore 

dn M 

dM — —b h (M\r)=l. (15) 
/ dM p m 



The | Jenkins et alJ ( 120011) mass function and Tink er et al.l 
( 120051) halo bias factors used here are fits to simulations over 
a finite range of halo masses, and they do not satisfy either 
of these constraints. To impose the constraint J 1 5I > explicitly, 
we break the second integral of equation il l\ at a halo mass 
Mbi-k = 10 8 h~ l M c .) and evaluate it as 



dM 2 



'0 
dn M 2 
dM 2 p m 

' dn M 2 
dM 2 — 

M brk dM 2 p m 



b h (M 2 \r) y m (k,M 2 )p no (x\Mi ,M 2 ) 



bh(M 2 \r) y m (k,M 2 )Pno(x\M 1 ,M 2 ) 



1- 



</M'-^7^h(M» 



MbA 



dM 1 



(16) 



For the term in brackets on the right hand side, we make the 
(good) approximation that p no (x) = y(k,M) = 1 for halos with 
M < Mbit at all radii of interest for our calculation, then ap- 
ply equation J 151 . We find that this procedure is necessary to 



obtain accurate results. The same problem does not arise for 
integrals involving (N)m because the mean occupation itself 
goes to zero at low halo masses. 

We must make one further adjustment to the analytic model 
before testing it against the populated A^-body halos described 
in § These halos are identified by the FOF algorithm, 
which, roughly speaking, selects particles within an isoden- 
sity surface. The mean overdensity A v i r within this surface 
depends on the halo profiles. To compute the effective value 
of A vrr for our simulations, we calculate the mean density 
within spheres centered on the most-bound particles of the 
FOF halos that enclose the halo's FOF mass. The canoni- 
cal value of A v ; r = 200 is accurate for our central model with 
Q m = 0.3 and a% =0.8. However, the A v j r values for other 
models, listed in Tabled deviate by up to 20%. The trend is 
as expected: halos in low Q m models are more concentrated 
because they form earlier, and they have higher values of A v ; r . 
However, the Q m -dependence of FOF halos selected with con- 
stant linking parameter is much weaker than the variation of 
virial densities predicted by the spherical collapse model (e.g., 
Brvan & Norman 1998). When calculating virial radii as a 
function of halo mass, we use the A v ; r values in Tabled 

The attached bottom frame in Figure |4^ shows the frac- 
tional difference between the analytic model and A^-body re- 
sults for the galaxy-galaxy correlation function, (Analytic ~ 
<£w-body) /6v-body, for the five cosmological/HOD models listed 
in Table Q Error bars represent the statistical uncertainty 
in the mean value of ^N-body, computed from the dispersion 
among the five independent simulations; for clarity, we show 
these only for the central model with £!,„ =0.3. The differ- 
ences between the analytic and numerical results are usually 
less than 5% at all r > 0.1 /r'Mpc. The one deviation that 
is clearly statistically significant is the rapid turn-up of the 
residuals at r ~ 0.1 /r'Mpc, which reflects the smoothing 
effect of the simulation's gravitational force softening. The 
marginally significant, ~ 5% discrepancy at r ~ 0.8 /r'Mpc 
suggests that our ellipsoidal exclusion correction still under- 
estimates the number of close halo pairs, even though it allows 
halos to be separated by less than the sum of their virial radii. 
Without this correction, the d eviation would fall off the bot- 
tom of the plot (see Fig. 10 in lTinker et all2005l) . The ~ 8% 
divergence of the model s at r ~ 8 /r'Mpc results from the 
deviations of the lSmith et alJ ll2003l) non-linear matter power 
spectrum from our simulation results. This difference could 
be partly an artifact of our finite box size, but the systematic 
dependence on cosmological parameters and reconvergence 
of results at r = 20 /r'Mpc suggest tha t it is mostl y a result 
of slight non-universality of the S mith et all ll2003l) formula, 
though it is hard to reach a definitive conclusion because of 
the substantial statistical uncertainties at these scales. 

The bottom frames in Figure 03 and|4j; show equivalent 
fractional differences for £ gm (r) and AS(r). The residuals 
for £gm(r) are similar to those for £gg(r), with the rise at 
r ~ 0.1 /r'Mpc, again reflecting the force softening in the 
numerical simulation. Since AE(r) depends on £ gm at all sep- 
arations less than r (see eq.Q), and the deviation between 
analytic and A^-body £ gm grows at smaller separations, a com- 
parison between the pure analytic calculation and the mea- 
surement of AS(r) from the simulation shows a substantial 
offset at r < 1 /r'Mpc. However, this offset reflects the lim- 
ited resolution of the A^-body simulation, not the failure of the 
an alytic model. (Similar dev iations at small scale were found 
by Mandelbaum et al. 2005.) To remove this numerical arti- 
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fact from the comparison, we set AE for the analytic model 
equal to the A^-body value at r = 0.1 /r'Mpc, then use the an- 
alytic calculation of £ em (V) to obtain the surface density for 
r > 0.1 /r'Mpc. With this correction, the analytic model for 
A£(r) is accurate to 5% or better for all five cosmological 
models at all radii 0.1 /r'Mpc < r < 20 /r'Mpc. 

4.2. Dissection of Correlation Functions 

It is interesting to examine the separate contributions of 
central and satellite galaxies to the galaxy-galaxy and galaxy- 
matter correlation functions. Figure |5p shows the familiar 
decomposit ion of ^(r) into 1 -halo and 2-halo contributions 
(see B erlind & Weinbergll2002l for extensive discussion). We 
adopt the fiducial model with il m = 0.3, as = 0.8, and the cor- 
responding HOD parameters listed in Table ^ and we use 
the A^-body simulation measurements. The 1-halo term dom- 
inates at small scales, but it drops rapidly towards larger r as 
halos with large virial diameters become increasingly rare. At 
large scales, the 2-halo term is a linearly biased version of the 
matter correlation function, but at small scales it turns over 
and drops because of halo exclusion. We plot 1 + £(r) rather 
than £(r) itself because pair counts are additive, so individual 
contributions sum to give the total 1 +£(r); this consideration 
becomes especially important for the central-satellite decom- 
positions shown in subsequent panels. The transition between 
1-halo and 2-halo dominance occurs at roughly the virial di- 
ameter of halos, where M* is the characteristic scale of the 
halo mass function. 

Figure [5J) separates 1 + £ gg (r) into central and satellite 
galaxy contributions. Central-satellite pairs dominate £ gg (r) 
at r < 0.4 /r'Mpc, and central-central pairs at r > 1 /r'Mpc, 
while satellite-satellite pairs dominate by a small factor in 
the intermediate regime. Figure |5J; and |5Jl show the sepa- 
rate central and satellite contributions to the 1-halo and 2- 
halo terms, revealing the origin of the behavior in Figure [5J). 
Central-satellite pairs dominate £' h (r) at small scales because 
in this regime most pairs come from the most common ha- 
los that are large enough to host a galaxy pair. These ha- 
los have (N)m < 3, and they therefore have more central- 
satellite galaxy pairs than pairs that involve only satellites. 
Conversely, satellite-satellite pairs dominate £ lh (r) at large 
scales because the halos with virial diameters large enough 
to host pairs at these separations have (N)m > 3. Finally, 
pairs involving at least one central galaxy dominate the 2- 
halo term by a large factor at all separations, because only 
halo pairs in which both halos have (N)m > 2 contribute, on 
average, more satellite-satellite pairs than central-central or 
central-satellite pairs. Such halos are much less common than 
those with 1 < (N)m < 2. Thus, satellite-satellite pairs make 
a major contribution to the total £ eg (r) only over the radial 
range in which 1-halo contributions from high mass halos are 
dominant. 

Figure |6] shows the equivalent dissection of £ em (r). The 
overall behavior is very similar to that seen in Figure[5] The 2- 
halo term extends to somewhat smaller scales because halos 
with mass near M m ; n can still form galaxy-matter pairs with 
lower mass halos that have smaller virial radii. The satel- 
lite contribution to the 2-halo term is analogous to the sum of 
central-satellite and satellite-satellite pairs for £gg( r ) because 
there are no "central" dark matter particles. (The normaliza- 
tion is lower than in Figure |2]because we now calculate the 
expected number of galaxy-matter pairs using all galaxies in- 
stead of satellites alone.) However, it is still a factor of 3-5 



below the central 2-halo term at nearly all separations. From 
Figure[6] we can understand why the individual halos of satel- 
lite galaxies appear to have so little impact on the SPH results 
discussed in §|2] Satellite galaxies dominate the 1-halo term 
°f £gm(r) only beyond r ~ 0.25 /r'Mpc. However, the indi- 
vidual halos of satellites orbiting in larger groups are usually 
tidally truncated well inside this radius, on scales where the 
signal is swamped by the contribution from central galaxies. 
Note that the satellite fraction of galaxy samples is typically 
less than 30%, and hence the lensing signal of satellite galax- 
ies is smaller by an order of magnitude than that of central 
galaxies. The dark halos of satellites in gro ups and clusters 
can b e measured by galaxy-galaxy lensing (Nataraian et al. 
120021) . but only by first identifying satellites and measuring 
A£(r) for them specifically. 

5. FROM GALAXY-GALAXY LENSING TO COSMOLOGICAL 
PARAMETERS 

Having established the accuracy of the analytic model, we 
can now use it to investigate the dependence of AE(r) on 
fi m and a%. We consider a well-defined sample of galax- 
ies, choose HOD parameters for each (J7 m , ag) combination 
by fitting the mean space density and projected correlation 
function of this sample, then predict A£(r). In this sec- 
tion, we focus on the sample of SDSS galaxies with M r < 
—21, with the projected correlation function, error covariance 
matrix, and mean space density n g = 1.17 x 10~ 3 /i 3 Mpc~ 3 
taken from Zehavi et al. (2005b). We also present predic- 
tions for a fainter lum inosity threshold, M r < -20, again using 
Zehav iet alJ i2005bh , s observational constraints. 13 At large 
scales, we expect to recover the linear theory, linear bias re- 
sult, AE oc Sl m /b oc <Tgf2 m , but we can extend the predictions 
to intermediate and small scales using the full analytic model. 

We make two significant changes in our application of the 
analyt ic model. First, we use a cmbfast transfer function 
(Seliak & Zaldarriaga 1996), comp uted for VL m = 0.3, h = 0.7, 
£lb = 0.04, in place of the Efstathiou et al. ( 1992) parametriza- 
tion adopted in our A^-body simulations. This change to the 
transfer function has little effect on the HOD parameters in- 
ferred by fitting w p (r p ), but it has a noticeable effect on the 
X 2 values of these fits, and it affects the AE(r) predictions 
themselves. Note that we do not change the transfer function 
when changing Cl m from our fiducial value of 0.3; because the 
power spectrum shape is empirically well constrained, we as- 
sume that any effect of changing ft m will be compensated by 
adjusting h, ilb, or the inflationary index n (which we set to 
one). Second, we define halo virial radii assuming A v j r = 200 
for all O m , instead of the varying A v i r values listed in Ta- 
bleland used in §|4] This change simply amounts to a slight 
change in the halo definition; to identify these halos in A^-body 
simulations, one woul d need to adjust the F OF linking length 
slightly with tt m . The Bullock et all (120011) concentration pa- 
rameters are defined for different (f2 m -dependent) A v i r values, 
but we rescale them to our A v i r definition. We still adopt the 
Jenki nset alJ yPOl) halo mass function for all cosmological 
models, with no rescaling. 

FigureQillustrates the results for a sequence of models with 
(78 ranging from 0.6 to 1.0, all fo r VL m = 0.3. T he HOD param- 
eter values required to match the Zehavi et al. ( 2005b) w p (r p ) 
measurements are listed in Table [2] and the mean occupation 
functions (N)m are shown in panelEfr. F° r lower a%, matching 

13 Specifically, we use the j Zehavi et al. 1 2005b) measurements for the 
M r < -20 sample with limiting redshift z — 0.06. 
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FIG. 5. — Dissection of the galaxy-galaxy correlations for the fiducial model, (a) Contributions of one-halo (dashed) and two-halo (dotted) galaxy pairs to the 
full correlation function (solid), (b) Contributions of central-satellite (dashed), satellite-satellite (dotted), and central-central (dot-dashed) galaxy pairs. Panels 
(c) and (d) show the central/satellite decompositions of the one- and two-halo terms individually. 



the observed clustering requires a larger fraction of galaxies in 
more massive, more biased halos, hence higher values of the 
(^Vsat)M slope a sat . The resulting galaxy correlation functions 
are very similar for all five values of a%, as shown in Figure|7Jl. 
Figure^p plots (M/p m ) dn/dlnM, the fraction of mass con- 
tained in a (natural) logarithmic bin centered at mass M. For 
(78 = 0.6, this function peaks near M ~ 10° /t'Mq, while for 
erg = 1-0 it peaks near M ~ 10 14 h~ l M Q . FigureQ; plots the 
same function multiplied by (N)m, a product that is propor- 
tional to the number of 1 -halo galaxy-matter pairs that arise in 
halos of mass M. This function peaks at M ~ 3 x 10 14 /i _1 M 
for cr 8 = 0.6 and M ~ 10 15 h~ l M Q for erg = 1.0. The trend of 
a sa t with us partly compensates the trend of the mass distribu- 
tion in Figure^, reducing the order-of-magnitude shift in the 
peak location to a factor of three. While high mass halos near 
the peak contribute a substantial fraction of all galaxy-matter 
pairs, these pairs are spread over a larger projected area, so 
the contribution in a given r ~ r + dr bin is multiplied by an 
additional factor that scales roughly as R~f r ~ M" 2 / 3 . For all 
values of a%, the fraction of 1-halo galaxy-matter pairs is tiny 
forM>5 x 10 15 h~ l M Q . 



TABLE 2 

HOD Parameters for Different <t 8 Models 



Model 




OS 


M ram (/i-'M ) 


M 1 (r 1 M ) 


asat 


1 


0.3 


0.6 


4.04 x 10 12 


6.28 x 10 13 


1.52 


2 


0.3 


0.7 


4.46 x 10 12 


7.98 x 10 13 


1.40 


3 


0.3 


0.8 


4.71 x 10 12 


9.58 x 10 13 


1.31 


4 


0.3 


0.9 


4.85 x 10 12 


1.11 x 10 14 


1.25 


5 


0.3 


1.0 


4.95 x 10 12 


1.23 x 10 14 


1.19 



NOTE. — The HOD parameters are chosen to reproduce 
the same clustering of the SDSS galaxy sample of M r < — 21 
and to match the number density of galaxies n g = 1.17 X 

10 -3 (/r'Mpcr 3 . 



Figures^ and[7J show the galaxy-matter correlation func- 
tions and excess surface density profiles, respectively, for this 
model sequence. At large scales, £, gm (r) and AE(r) increase 
with <78 oc 1 /b as expected from linear theory. A similar in- 
crease appears on small scales because of the larger fraction of 
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FIG. 6. — Dissection of the galaxy-matter correlations for the fiducial model, in the same format as Fig. |5] In panel (b)-(d), dashed and dotted lines show 
galaxy-matter pairs involving a central galaxy or a satellite galaxy, respectively. 



galaxy-matter pairs in more massive halos. In fact, the shapes 
of £ gm (r) and A£(r) appear remarkably constant over the full 
range 0.1 /r'Mpc < r < 20 /r'Mpc, a point we quantify be- 
low. 

Figure [8] shows A£(r) for three model sequences with dif- 
ferent variations. In Figure^, we consider the same sequence 
of increasing erg, fixed Q m shown in Figure Q but we al- 
ways keep halo concentrations fixed at the values predicted 
for erg = 0.8. We adjust HOD parameters slightly from the 
values listed in Table [2] to obtain the minimum-x 2 fit to the 
projected correlation function with the new halo concentra- 
tions. In the large r, 2-halo regime, AE(r) is nearly identical 
to that shown in Figure [7J . However, fixing the concentra- 
tions to those of the central model has an important effect 
at small scales, causing the AE(r) curves to converge. The 
constancy of shape in Figure is thus partly a consequence 
of the changes in halo concentrations in different erg models; 
higher a% leads to earlier halo collapse and higher concentra- 
tion, boosting AE(r). 

Figure [8p shows a sequence with fixed erg = 0.8 and ft,,, 
varying from 0.2 to 0.4 in steps of 0.05. For this sequence, 
the HOD parameters M m [ n and Mi scale in proportion to Q m , 



though we again make slight adjustments to fit w p (r p ). If 
halo concentrations were independent of fi m , and we did not 
make those small adjustments, then £gg( r ) and £gm(r) would 
be identical for all five models after this mass rescaling, since 
the halo mass function and ha lo bias factor s are functions of 
M/M* and M* oc fl m ( see lZheng et al.l2002l for further discus- 
sion). In this case, AS(r) would have a constant shape and an 
amplitude proportional to f2,„. Figure HJ) shows roughly this 
behavior, but the trend of higher concentration for lower Q. m 
produces a weak convergence of models at small r. 

We have so far assumed that satellite galaxies trace the dark 
matter in halos, with the same NFW radial profile. We now 
relax this assumption and allow the satellite profiles to have 
a lower concentration param eter, as suggested by som e nu- 
merical studies (see §[2] and Nagai & Kravtsov 2005). Fig- 
ure|SJ; compares models with il m = 0.3, erg = 0.8, and satellite 
concentration parameters c ga i = a c c±a with a c = 0.3, 0.5, 0.7, 
and 1.0. We again adjust HOD parameters to fit w p (r p ) af- 
ter changing galaxy concentrations. These adjustments partly 
compensate for the changes in galaxy concentration, so the 
effect of a radial profile change is somewhat smaller here 
than in § |2 (Fig. 0, where we kept other HOD parame- 
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FIG. 7. — Clustering contributions and clustering signals for the model sequence with fixed Q m and varying <jg. (a) Mean halo occupation functions, determined 
by fitting the w p (r p ) data, with <r 8 increasing from top to bottom (see panel d legend), (b) Fraction of matter per logarithmic bin of halo mass, (c) Fraction of 
galaxy-matter pairs per logarithmic bin of halo mass. Note that pairs in higher mass halos are spread over a wider range of separations, diluting the contribution 
to any given r — > r + dr bin. Panel (d), (e), and (f) show the galaxy-galaxy correlation function, galaxy-matter correlation function, and excess surface density, 
respectively. 
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FIG. 8. — Excess surface density profiles for other model sequence. We consider (a) the same model sequence as in Fig.Q but with halo concentrations held 
fixed at the values for Q m = 0.3 and erg = 0.8, (b) a sequence of models with fixed <rg = 0.8 and Q m ranging from 0.2 (dotted) to 0.4 (long dashed) in steps of 
0.05, and (c) models with Q m = 0.3 and a% = 0.8 in which galaxy profile concentrations satisfy Cgai = a c Cd m , with a c = 0.3 (dotted), 0.5 (long-dashed), 0.7 
(dot-dashed), and 1.0 (solid). The inset panel shows fractional deviations from the a c = 1.0 model. 



ters fixed. For a c > 0.7, the impact on AS(r) is under 3% 
at all r. For a c = 0.3, the effect rises to 7% at the scale 
r ~ 0.5- 1 /i _1 Mpc where satellite galaxies make a domi- 
nant contribution to £ gm (r) (see Fig. |(J>). Small differences 
between galaxy and dark matter concentrations can thus be 
safely neglected, but large differences can have a small but 
measurable impact at r < 2 /i _1 Mpc. 

Observational studies of galaxy-galaxy lensing often exam- 
ine the ratio £. Pm {r) I £ vv (r), which is equal to r„m/b in th e lin- 
ear bi as model (e.g.. Hoe kstra et aH 2001 . 2002t lSheldon et all 
2004). Figure [5] plots this ratio for the four model sequences 
shown in Figures and [8] Figure^ shows the fixed-f2 m , in- 
creasing a% sequence of Figure Since HOD parameters 
in each model are adjusted to match the projected correla- 
tion, model differences are driven almost entirely by changes 
in £ gm (r). At large scales, the ratio £ gm (r) / £,gg( r ) is con ~ 
stant as predicted for linear bias, and t; gm (r) increases in 
proportion to a% oc 1 jb. Comparison to the bias factor de- 
fined by b 1 = Cgg( r )/Cmm('") implies a cross-correlation co- 
efficient r gm ~ 0.9 for all five models. However, in every 

case £gm( r )/£gg( r ) rises sharply at a scale r ~ 1 /r'Mpc near 
the 1-halo to 2-halo transition. For £ gg (r), this transition is 
fairly shar p, producing measura ble deviations from a power- 
law JBerhnd & Weinberg|l200l iZehavi etalJl200l . These 
deviations are smoothed out in £ gm (r) because the contribu- 
tion of halos below M m ; n allows the 2-halo term to overlap 
more with the 1-halo term (compare Figs.|5^ and[(J), and the 
ratio ^gmW/Cgg( r ) therefore shows a sharp feature reflecting 
the break in £ gg (r). For higher cr%, the 1-halo term extends 
to larger r, and the jump in £gm( r )/£gg( r ) is larger in ampli- 
tude but spread over a larger range of r. The break in £ gg (r) 
is generall y stronger for more strongly c lustered galaxy sam- 
ples jBerlind & Weinbergl2002HZehavi et alJ2005 b). and we 
expect similar sample dependence for the £, gm (r) /£ gg (r) jump. 

At small scales, £ gm (r)/£gg( r ) is again roughly flat, but at 
a level higher than the large scale ratio. As noted earlier, 
fixing halo concentrations causes the galaxy-matter correla- 



TABLE 3 
The Weak Lensing Signal 





M r < 


-21 




M r < 


-20 




logr 


logAiJ(r) 


a 


H 


logAlV) 


a 


P 


-1.029 


1.929 


0.65 


1.13 


1.582 


0.71 


0.84 


-0.952 


1.870 


0.68 


1.10 


1.513 


0.73 


0.83 


-0.875 


1.808 


0.71 


1.08 


1.446 


0.75 


0.83 


-0.798 


1.743 


0.73 


1.06 


1.378 


0.77 


0.83 


-0.721 


1.674 


0.76 


1.04 


1.310 


0.79 


0.85 


-0.643 


1.602 


0.79 


1.03 


1.244 


0.81 


0.85 


-0.567 


1.529 


0.80 


1.03 


1.181 


0.83 


0.87 


-0.489 


1.453 


0.83 


1.03 


1.116 


0.84 


0.91 


-0.412 


1.379 


0.85 


1.02 


1.050 


0.86 


0.96 


-0.335 


1.306 


0.87 


1.03 


0.982 


0.88 


1.00 


-0.258 


1.227 


0.89 


1.06 


0.922 


0.89 


1.08 


-0.181 


1.145 


0.91 


1.10 


0.860 


0.91 


1.13 


-0.104 


1.064 


0.92 


1.15 


0.806 


0.92 


1.16 


-0.027 


0.983 


0.93 


1.19 


0.753 


0.95 


1.19 


0.050 


0.891 


0.94 


1.24 


0.686 


0.96 


1.24 


0.127 


0.793 


0.95 


1.31 


0.610 


0.98 


1.31 


0.204 


0.697 


0.96 


1.37 


0.535 


0.99 


1.38 


0.281 


0.604 


0.98 


1.41 


0.456 


1.00 


1.44 


0.358 


0.511 


0.99 


1.45 


0.378 


1.01 


1.49 


0.435 


0.419 


1.00 


1.48 


0.297 


1.02 


1.53 


0.512 


0.330 


1.00 


1.49 


0.218 


1.03 


1.55 


0.589 


0.244 


1.01 


1.47 


0.137 


1.04 


1.55 


0.666 


0.162 


1.02 


1.43 


0.059 


1.04 


1.51 


0.743 


0.089 


1.03 


1.37 


-0.013 


1.05 


1.43 


0.820 


0.021 


1.04 


1.29 


-0.076 


1.05 


1.35 


0.897 


-0.039 


1.03 


1.21 


-0.131 


1.06 


1.27 


0.974 


-0.094 


1.05 


1.14 


-0.196 


1.05 


1.20 


1.051 


-0.145 


1.05 


1.09 


-0.242 


1.05 


1.15 


1.128 


-0.193 


1.05 


1.06 


-0.288 


1.05 


1.12 


1.205 


-0.241 


1.04 


1.05 


-0.333 


1.04 


1.12 


1.282 


-0.292 


1.04 


1.05 


-0.382 


1.04 


1.12 



NOTE. — The excess surface densities ASfid of the fiducial model (Q m = 
0.3, erg = 0.8) for SDSS g alaxy samples of M r < -21 and M, < -20. The 
projected radius r is in /r'Mpc, and AS is in hMQpc" 2 . 



tion function of different a% models to converge (Fig. [8^), so 
the c; gm (r)/c; gg (r) ratios also converge in this case (Fig.|9p). 
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FIG. 9. — The ratio £gm(r)/£g g (>), which is equal to r gm /b in the linear bias model. Panels (a)-(d) show the four model sequences illustrated in Figs. 171 andlSl 
Line types follow the same sequence as in those Figures, with erg increasing from 0.6 (dotted) to 1.0 (long-dashed) in panels (a) and (b), Q„, increasing from 
0.2 (dotted) to 0.4 (long-dashed) in (c), and a c increasing from 0.3 (dotted) to 1.0 (solid) in (d). 



Figures |5J; and |5Jl show that the effects of fi m or c ga i varia- 
tions are much smaller than those of erg variations. The ~ 5% 
model-to-model differences at small scales arise from their 
different halo concentrations, while the smaller differences at 
large scales reflect the slight changes in HOD parameters re- 
quired to match w„{r„). P resent observations (iHoekstra et al.l 
120021 ISheldon et aljl2004l) are consistent with £g m (r)/£ gg (r) 
that is approximately scale-independent, but the uncertainties 
are still fairly large, and testing for the feature predicted in 
Figure|9]will require more careful replication of observational 
procedures. 

Figure [9] shows that the linear bias expectation of constant 
Cgm('")/^gg('") holds accurately for r > 4 /r'Mpc but fails at 
the 20-50% level in the non-linear regime. We can also ask 
how well the linear bias prediction AS oc Vl m <j% describes the 
scaling of AE(r) with cosmological parameters. To answer 
this question, and to allow easy scaling of our predictions with 
cosmological parameters, we adopt the more general formula 

AS(r) _fn m \ a ^a & Y , m 
AS F i D (f) \0.3J V0.8/ 
and determine best-fit values of a and (3 at each separation 



r. Here AS F i D (r) is the excess surface density prediction of 
the fiducial model with Cl m = 0.3 and erg = 0.8, and we fit a 
and f3 using a full grid of models with erg = 0.6, 0.7, 0.8, 0.9, 
1.0 and Ct m = 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45. We assume 
Cgai = cam in all cases. 

Figures flOk and 1 10b plot the fitted values of a and (3, re- 
spectively, as a function of r. Results for M r < -20 galax- 
ies and M r < -21 galaxies are similar, though the underlying 
ASfidM i s different in the two cases. There are two notable 
departures from the linear bias values a = f3 = 1. At small 
scales, a falls below one, reflecting the weak convergence of 
AS(r) curves seen in Figure|8j5. This convergence in turn re- 
flects the il m -dependence of halo concentrations. At scales 
r ~ 2-5 /r'Mpc, (3 rises above unity, corresponding to the 
slight divergence of AE(r) curves at these scales in Figure[7F. 
FigureflOfc shows the rms and maximum fractional errors be- 
tween the ft,,, or erg dependences predicted by the full ana- 
lytic model and the scaling relation (I17t . calculated over our 
full model grid. The rms errors range from ~ 1% at large 
r to ^ 3% at intermediate r. The largest errors arise for the 
Q m = 0.45, erg = 1.0 model, and they are roughly twice the 
rms errors. Figure ITOh shows the result of adopting the linear 
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FIG. 10. — Parameters a(r) and (3(r) of the bias scaling relation AE(r) oc Sl°(j^ (eg. 1171 ) for samples matched to the SDSS w p (r p ) measurements of M r < —20 
(triangles) and M r <—2i (circles) galaxies. Panel (c) shows the rms and maximum fractional errors of this scaling relation, relative to the full analytic calculation, 
over a model grid with Q m varying from 0.15 to 0.45 and erg varying from og = 0.6 to 1.0. Panel (d) shows the same errors for linear bias scaling a = (3 = 1. In 
all cases the maximum error arises for Q m = 0.45, erg = 1.0. 



bias scaling a = (3 = 1 in equation d!7i . The linear bias pre- 
dictions are accurate at the ~ 1% (rms) to ~ 3% (maximum 
error) level for r > 10 /r'Mpc, but the deviations become sub- 
stantial at smaller r, with errors of ~ 5- 16% at r = 3 /r'Mpc 
and r = 0. 1 /i _1 Mpc for the Cl m = 0.45, cr 8 = 1 .0 model. The er- 
rors of the linear bias scaling are typically a factor ~ 1 .5-4.0 
larger than those using our fitted values of a(r) and (3(r). 

Table |3]lists AfidM' me values of AS(r) predicted by the 
analytic model for the M r < -20 and M r < -21 galaxy sam- 
ples assuming fi m = 0.3 and a% = 0.8. It also lists the a(r) 
and j3(r) functions shown in Figure Equation ( 1171 can 
be used to scale these predictions to other values of f2 m and 
<j%, and measurements of AE(r) for these galaxy samples 
could then be used to obtain constraints in the f2 m -(T8 plane. 
Our prediction of ASfidM is weakly dependent on the HOD 
parametrization that we adopt when fittin g w,,(r n )- If we adop t 
the alternative parametrization used by Zeh avi et alJ J2005b). 
where a sat is fixed to one but the (N sat )M cutoff is a fit param- 
eter, then the AS(r) predictions change by less than 5% for 
all r > 0.1 /i _1 Mpc. We have not yet explored more general 
parametrizations, but we expect that the AS(r) predictions 



would be robust at this level for all HOD models that fit the 
observed w p (r p ). 

Galaxy-galaxy lensing measurements are often made for 
flux-limited samples rather than absolute-magnitude limited 
samples to increase the signal-to-noise ratio, but such mea- 
surements are difficult to interpret quantitatively because they 
do not represent properties of a uniformly defined galaxy 
population. Because our predictions apply on the non-linear 
scales where the measurement precision is higher, and be- 
cause results from different radii can be combined, it should 
be possible to obtain precise constraints on crgQ m from abso- 
lute magnitude-limited samples, and to use different samples 
to check for consistency. Since the values of a and j3 vary 
with scale, it is possible in principle to break the degeneracy 
between Q m and erg. However, the deviations from linear scal- 
ing are not large, so while a%£l m should be well constrained, 
individual parameter constraints are likely to be imprecise and 
s ensitive to systemati c uncer tainties in the modeling. 

Mandelbaum et al. (2006) have recently presented SDSS 
measurements for narrow bins of luminosity and of stellar 
mass, which are well suited to their goal of constraining halo 
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virial masses and satellite fractions. For cosmological param- 
eter constraints, we think it is better to use luminosity or mass- 
threshold samples, which provide higher signal-to-noise ratio, 
and which are easier to model robustly because there is no 
upper mass cutoff on (N cen )M- 

6. SUMMARY 

We have developed an analytic model to predict A£(r) 
for specified cosmological and galaxy HOD parameters and 
tested its validity using SPH and A^-body simulations. We 
have used the analytic model to investigate the dependence 
of AE(r) on erg and Vl m when HOD parameters are chosen to 
reproduce the observational space density and projected cor- 
relation function w p (r p ) of the galaxy sample being measured. 
Our main findings are as follows: 

1 . In our SPH simulation, replacing the satellite galaxies of 
each halo with randomly selected dark matter particles has a 
10-20% effect on £ gg (r) and £ gm (r) at scales r ~ 0.5 /r'Mpc, 
and smaller impact at other scales. Most of this difference 
arises from the differing radial profiles of satellite galaxies 
and dark matter. If satellites are replaced in a way that pre- 
serves the radial profile but randomizes azimuthal positions, 
then changes to £ gg (V) and £ gm (/) are < 10% at all radii, and 
changes to A£(r) are < 5%. Dark matter subhalos around in- 
dividual satellites orbiting in larger halos are present, but they 
have negligible impact on the global AE(r). 

2. If we randomly reassign the galaxy occupation num- 
ber of each halo with M < 4.6 x 10 13 h M® to another halo 
of nearly equal mass, then changes to £gg( r X £gm( r )> an d 
AS(r) are < 2% at all r < 5 /r'Mpc. This result implies that 
any environmental dependence of the halo occupation func- 
tion P(N\M) at fixed halo mass has minimal impact on these 
statistics for our simulated galaxy sample, which is defined 
by a baryonic mass threshold. For our largest scale point at 
12 /r'Mpc, we find an effect of 10% on £ gg (r), 5% on £ gm (r), 
and 2% on A£(r), but the statistical uncertainties of our esti- 
mate are of comparable magnitude at this scale, so larger sim- 
ulation volumes are needed to definitively establish the impact 
of any environmental dependence on the large scale bias fac- 
tor. Taken together, results 1 and 2 show that the £gg( r ) and 
AS(r) predictions of a full hydrodynamic simulation can be 
reproduced to 5% or better (usually much better) by populat- 
ing the halos of a pure A^-body simulation with the correct 
P(N\M), provided that satellite galaxy populations have the 
correct radial profiles. 

3. Our analyt ic model for A£ (r) is based on the me thods 
introduced bv lSeliakl (120001) and lGuzik & Seliakl (120011) . but 
it incorporates the scale-dependent halo bias and ellip soidal 
halo exclusion corrections introduced bv lZhenglll200l and 
Tinker et al. (2005) for £ gg (V) calculations. We have tested the 
analytic model against numerical results from a grid of pop- 
ulated A^-body simulations, which span the parameter range 
cr 8 = 0.6-0.95 and O m = 0.1-0.63, with HOD parameters 
chosen to match the space density and projected correlation 
function of SDSS galaxy with M- < -20 ( Ze havi et alJ2 005b). 
The analytic model reproduces the numerical results to 5% or 
better over the range of 0.1 /r'Mpc < r < 20 /r'Mpc. The 
residuals are consistent with the statistical errors of the numer- 
ical calculations, except for the innermost bin, where gravita- 
tional force softening in the A^-body simulations artificially 
suppresses correlations. 

4. For the M r < -20 HOD parameters, pairs involving at 
least one central galaxy dominate the galaxy-galaxy correla- 
tion function at r < 0.4 /r'Mpc and r > 1 /r'Mpc. In the 



range 0.4 /r'Mpc < r < 1 /r'Mpc, satellite-satellite pairs 
in large halos make the dominant contribution to £ gg (r). In 
similar fashion, central galaxies dominate the galaxy-matter 
correlation function at small and large separations, while 
satellite galaxies dominate in the range 0.25 /r'Mpc < r < 
1.5 /r'Mpc. The halos of individual satellites make negli- 
gible contribution to AS(r) because they are usually tidally 
truncated below the scales at which the satellite contribution 
itself is important. 

5. For samples with HOD par ameters chosen t o match 
the M r < -21 SDSS sample of lZehavi etafl l l2005bl) . the ra- 
tio £gm(r)/£gg(V) is constant at r > 4 /r'Mpc, as predicted 
by the linear bias model, but it jumps by 20-50% at scales 
r ~ 1 /r'Mpc near the transition from the 1-halo to 2-halo 
clustering regime, before settling to a new, higher value at 
small scales. The magnitude of the jump depends on cr%, 
and it is likely to depend on the galaxy sample as well, be- 
ing stronger for more highly clustered populations. In linear 
bias terms, the large scale values of £, gm (r) / £, sg (r) correspond 
to a galaxy-matter cross-correlation coefficient r gm ~ 0.9, if 

we define the bias factor b = [£ gg (r)/£mm (/■)]■ 

6. We fit the dependence of AS(r) on cosmological param- 
eters with a scaling formula AE(r) cx fi^Og , where a and (3 
are slowly varying functions of r. This scaling describes the 
results of our full analytic model with rms error < 3% over 
the parameter range Q m =0.15- 0.45, erg = 0.6 - 1 .0. At large 
scales, a and j3 approach the linear bias values a = (3 = 1 . 
However, forcing a = (3 = 1 at all scales leads to errors that 
are larger by factors of 1.5-4.0, relative to the scaling for- 
mula d!7l > with our fitted values of a and (3. 

Table[3]lists our predicted values of AS(r) for the M r < —20 
and M r < -21 SDSS samples, assuming our fiducial cosmo- 
logical model with fi m = 0.3 and erg = 0.8. Equation Jl 71 al- 
lows scaling of these results to other values of erg and f2 m , and 
measurements of AE(r) for these samples can be combined 
with these predictions to obtain cosmological constraints, 
which will be tightest on a parameter combination that is ap- 
proximately <7gf2 m . Given the grow t h of th e S PSS since the 
samples analyzed by ISheldon et alJ (120041) and IZehavi et alJ 
(2005b), it is probably preferable to extract w p (r p ) and AS(r) 
estimates for matched galaxy samples from the latest data 
sets. A full analysis should also investigate the effects 
of adding greater flexibility to the HOD parametrization it- 
self, u sing, e.g., the 5-parameter formulation of IZheng et alJ 
(2005). We have tested the effect of changing to a different 3- 
parameter description (see §[5jl, and we find changes of < 5% 
in the AS(r) predictions. We suspect that these predictions 
would remain similar for any choice of HOD parameters that 
reproduces the observed w p (r p ). Our SPH and A^-body tests 
indicate that the analytic model predictions should be accu- 
rate to 5% or better given our HOD parametrization, though 
assumptions about galaxy profile concentrations have signifi- 
cant effect at r < 2 /r'Mpc. This level of accuracy is adequate 
given the statistical errors expected for current samples, but 
refinement and testing on large simulations will be needed to 
take full advantage of future analyses of even larger, deeper 
surveys. Precise determinations of £l m and erg can play an im- 
portant role in testing theories of dark energy and models of 
inflation, making galaxy-galaxy lensing an essential element 
of observational cosmology. 
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